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J/ '4) production has been measured in proton-proton collisions at -^i = 200 GeV over a wide 
rapidity and transverse momentum range by the PHENIX experiment at RHIC. Distributions of 
the rapidity and transverse momentum, along with measurements of the mean transverse momentum 
and total production cross section are presented and compared to available theoretical calculations. 
The total J/tp cross section is 3.99 ± 0.61(stat) ± 0.58(sys) ± 0.40(abs) jih. The mean transverse 
momentum is 1.80 ± 0.23(stat) ± 0.16(sys) GeV/c. 

PACS numbers: 13.85.Ni, 13.20.Fc, 14.40.Gx, 25.75.Dw 



Understanding J ftp production mechanisms requires 
data over a large range of collision energies and with 
broad coverage in rapidity and transverse momentum 
(pt)- Existing data at lower energies from fixed tar- 



get hadron experiments yield total cross sections and 
mean pT ({pt)) values in the energy range y/s = 7 - 
38.8 GeV Limited kinematic coverage in collider ex- 
periments ^ |3j 01 has so far meant that total cross sec- 
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tions and mean pT values could not be measured. The 
systematic study of J/ ip production at Relativistic Heavy 
Ion Collider (RHIC) energies with wide pt and rapid- 
ity coverage should therefore provide crucial tests of J/ip 
production models. In addition, the RHIC proton-proton 
results provide a baseline for studying cold and hot nu- 
clear matter in proton-nucleus and nucleus-nucleus colli- 
sions using Jf-ip yields as a probe. 

Intense theoretical interest in the J/tp production 
mechanism was stimulated when the Color Singlet Model 
(CSM) was found Q to dramatically underpredict the 
high Pt CDF prompt J/ip and ip{2S) cross sections 0. 
Attention turned toward models in which color octet 
cc states can also contribute to the J/^ yield. The 
Color Octet Model (COM_K which is based on the Non- 
Relativistic QCD model pj, has been successful in re- 
producing the high px CDF prompt J/i/j cross sections, 
as has the more phenomenological Color Evaporation 
Model 0. 

In this paper we report results of the first measure- 
ments oipp J/ip+X at RHIC, made at = 200 GeV 
by the PHENIX experiment. The data yield the first to- 
tal cross sections for J/ ip production beyond fixed target 
energies, and the first measurement of (jjt) beyond ^/s 
= 63 GeV. They will constrain models in the lower px 
region where gluon fusion is expected to dominate (at px 
beyond about 5 GeV/c, the direct J/ip production cross 
section is expected to be dominated by fragmentation of 
high Pt gluons |3]). 

The PHENIX experiment Q detects electrons in the 
pseudo-rapidity range \r]\ < 0.35 in two central spectrom- 
eter arms covering = 90°, and forward rapidity muons 
in two muon arms covering Acj) = 360°. Only one muon 
arm, covering 1.2 < < 2.2, was operational for this 
data set. Electrons are identified by matching charged 
particle tracks to energy deposits in the Electromagnetic 
Calorimeter (EMC) and to rings in the Ring Imaging 
Cerenkov Detector (RICH), which has a threshold of 4.7 
GeV/c for pions. Muons are identified by finding deeply 
penetrating roads in the Muon Identifier (MuID) and 
matching them to tracks in the Muon Tracker (MuTr). 

The data were recorded during the 2001/2002 run at 
= 200 GeV. After quality assurance and vertex cuts 
(± 35 cm for ee and ± 38 cm for up), 67 nb~^ were used 
for the J/tp ^ A'^M^ analysis, and 82 nb~^ for J/tp ^ 

e~ . The minimum bias interaction trigger required 
at least one hit on each side of the interaction vertex in 
the Beam-Beam counter (BBC). Minimum bias trigger 
rates varied from 5 to 30 kHz. Events containing J/tp 
decays were selected using level- 1 triggers in coincidence 
with the minimum bias interaction trigger. The J/tp — > 

e~ trigger required a minimum energy deposit of either 
0.75 GeV in a 2 X 2 tile of EMC towers or 2.1 GeV 
in a 4 X 4 tile. The J/tp ^ /^^/-^^ trigger required at 
least two deeply penetrating roads in separate azimuthal 
quadrants of the MuID (9j. 
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FIG. 1: The invariant mass spectra for dielectron and dimuon 
pairs. Unlike-sign pairs are shown as solid lines, the sum of 
like-sign pairs as dashed lines. 



J/tp yields in the central arms were obtained by re- 
constructing electron-positron pairs. Electron candi- 
dates were charged particle tracks that were associated 
with a RICH ring (> 2 hit phototubes) and an EMC 
hit (±4(7 position association cut), and which satisfied 
0.5 < E/p < 1.5, where E is the EMC cluster energy 
and p is the reconstructed track momentum. A 5 GeV/c 
upper limit on electron momentum prevented charged pi- 
ons from firing the RICH. 

J/tp yields in the muon arm were obtained by recon- 
structing IJ,~^IJ,^ pairs. Muon tracks were reconstructed 
by finding a track seed in the MuID and matching it to 
clusters of hits in each of the three MuTr stations. The 
momentum was determined by fitting, with a correction 
for energy loss, the MuID and MuTr hit positions and the 
BBC vertex position. Each track was required to pass 
5(7 cuts on the from the track fit and on the radial 
distance of the fitted track from the measured z-vertex 
position. 

Unlike-sign pairs and, for background estimation, like- 
sign pairs satisfying the above conditions were combined 
to form invariant mass spectra. Simulations show that 
the acceptance for like-sign and unlike-sign pairs is the 
same to within a few percent for invariant masses above 
1 GeV/c^ for electrons. In Fig. ^ unlike-sign and like- 
sign invariant mass spectra from the entire pp data set are 
shown together. For electrons, the net yield inside 2.8 - 
3.4 GeV/c^ is 46, for muons inside 2.71 - 3.67 GeV/c^ it is 
65. For electrons, the peak width is 110 MeV/c^ and the 
centroid agrees well with the PDG value lo']. For muons, 
the width is 160 MeV/c^. The muon peak centroid is 
higher than the PDG value by about 3%, consistent with 
the uncertainty in the muon magnetic field calibration. 

The J/tp cross sections were determined from the mea- 
sured yields using 
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where Njj^ is the measured J/V' yield, / Cdt is the inte- 
grated luminosity measured by the minimum bias trigger, 
Bii is the branching fraction for the J/tj) to either e+e~ 
or pairs (PDG average value 5.9% 0), e^ms is the 

minimum bias trigger efficiency for an event containing 
a J/V', e-ivii is the level-1 trigger efficiency for detecting 
a J/V', and Acrec is the acceptance times reconstruction 
efficiency for a J /"ip- 

The integrated luminosity can be written as / Cdt = 
Nmb/o'bbc, where Nmb is the number of minimum bias 
triggers and ubbc is the minimum bias trigger cross sec- 
tion. Using a van der Meer scan measurement, ubbc 
was determined to be 21.8 ± 2.1 mb 0|. We have 
estimated euas in two ways. First, the minimum bias 
trigger efficiency for J /ip events from a simulation study 
using Pythia [i3| (with the GRV94NLO parton dis- 
tribution functions (PDFs)) was 0.74, with no pr de- 
pendence. Good agreement is observed in the dNch/drj 
distribution between Pythia simulations and measure- 
ments [l^ for events involving high pT tt'^ production. 
Second, we measured our minimum bias trigger effi- 
ciency for high pt tt'^ production using events recorded 
with a high pr EMC trigger. The efficiency of 0.75 ± 3% 
is constant within uncertainties over the measured range 
of 1.5 < Pt < 9 GeV/c. We chose to use our measured 
trigger efficiency from high pT 7r° events when calculating 
the J/ tp cross section, assuming that the trigger efficiency 
is the same for both processes. 

For the electron analysis, A Srec was determined as a 
function of px using a full GEANT simulation of single 
J/ip events with flat distributions in dN/dy {\y\ < 0.6), 
Pt {pt < 6 GeV/c) and coUision vertex {\z\ < 35 cm). 
The GEANT simulations were tuned to match real detec- 
tor responses for single electrons. The reconstruction effi- 
ciency calculations used a typical dead channel map. An 
average correction for run-to-run variations in detector 
active area was determined from single electron yields. 
An estimate of the systematic uncertainty in A Crec due 
to z vertex dependence of the acceptance, momentum 
resolution effects, the pair mass cut, and electron identifi- 
cation cuts is given in Tabled along with the uncertainty 
in the yield due to Drell-Yan and correlated charm decay 
contributions. The efficiency eivn of the level-1 J/V' trig- 
gers in the central arms was determined as a function of 
Pt by using a software trigger-emulator to analyze simu- 
lated single J/ip events. The results are shown in TableQ] 
The trigger emulator was tuned by analyzing simulated 
single electrons and comparing with the real single elec- 
tron trigger efficiency. 

For the muon arm, Aerec^ivii was determined as 
an average within each rapidity and pt bin, using a 
full GEANT simulation with J/ip events generated by 
Pythia (with GRV94LO PDFs). The Pythia J/ip ra- 
pidity and Pt distributions are very similar to those of the 
real data, so that bin averaging effects should be approx- 
imately accounted for by this procedure. The simulated 



TABLE L Table of quantities and their systematic error esti- 
mates. Ranges are given for pt dependent quantities. For the 
fj,^ case the values of Aerec and eivii axe combined. The ab- 
solute cross section normalization uncertainty from etias and 
/ jCdt is kept separate and is labeled (abs). 



Quantity 



Yield ± 5% 

Aer^c 0.026-0.010 ± 13% 

eiM (2x2) 0.87-0.90 ± 5% 
(4x4) 0.30-0.74 ± 36% 

euas 0.75 ± 3% 

/ jCdt 82 nh-^ ± 9.6% 

Total ± 15%(sys) ± 10%(abs) ± 14%(sys) ±10%(abs) 



± 5% 
0.038 - 0.017 ± 13% 



0.75 ± 3% 
67 nb-i ± 9.6% 
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FIG. 2: The J/i/i pr distributions for the dielectron 
and dimuon measurements, with statistical uncertainties. 
The solid line is a phenomenological fit of the form 
1/{2tvpt) da/dpT = A{1 + {pr/Bfy^. The dashed line is 
an exponential fit. The CSM (dot-dashed) and COM (dotted) 
calculations are from 15]. 



events were reconstructed using the same reconstruction 
software and cuts as for the real data, assuming nomi- 
nal detector efficiencies and typical realistic dead channel 
and dead high-voltage maps. Each event had to pass the 
simulated dimuon trigger. The systematic error includes 
discrepancies between Monte Carlo and real detector re- 
sponse, run to run variations in the detector state, and 
uncertainties in the Pythia distributions. The results, 
integrated over the rapidity range of the muon arm, are 
shown in Tabled along with an estimate of the system- 
atic error on the yield due to the background subtraction 
technique. In both the electron and muon cases the Jf-ip 
polarization was assumed to be zero, since existing J/^ 
polarization measurements are consistent with zero at 
low Pt 14] . The effect of the unknown J/^ polarization 
has not been included in the systematic error. 

The Pt distributions for J/ip e+e" and J/ip ^ 
1^.^1-1.^ are shown in Fig.[51 with predictions |l5l ] from the 
COM. Predictions of the CSM, which greatly underpre- 
dicts the cross sections, are also shown. These predictions 
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are limited to > 2 GeV/c because parton intrinsic 
transverse momentum (kx) broadening is not accounted 
for properly in the calculation. The COM calculations do 
not include fragmentation contributions, which become 
important at around 5 GeV/c j3- Calculations cover- 
ing all pt and including fragmentation contributions are 
needed. The solid lines are a phenomenological fit of a 
form that has been shown to fit J/ip data well at fixed 
target energies The dashed line is an exponential 

fit. The phenomenological fits yield (pt) values of 1.85 
± 0.46(stat) ± 0.16(sys) GeV/c (central arm) and 1.78 ± 
0.27(stat) ± 0.16(sys) GeV/c (muon arm), with a com- 
bined value of 1.80 ± 0.23(stat) ± 0.16(sys) GeV/c. The 
systematic uncertainties were estimated from the spread 
in (pt) from a weighted mean of the binned data, the 
phenomenological fit, and the exponential fit. An ad- 
ditional 3% was assigned to the muon (px) due to the 
uncertainty in momentum scale. 

The J/tp rapidity distribution obtained by combining 
the dielectron and dimuon measurements is shown in 
Fig. 121 with the muon arm data divided into two ra- 
pidity bins. The COM curves are theoretical shape pre- 
dictions using the same models as are discussed in 
connection with Fig. except that they are normal- 
ized to our data to make the shape comparison clearer. 
Since gluon fusion is the dominant process in all of the 
models, the rapidity shape depends mostly on the gluon 
distribution function and is not very sensitive to the pro- 
duction model. Most of the available PDFs are consistent 
with the data, and improved statistical precision will be 
needed to constrain them. A Pythia calculation that 
reproduces the shape of our data best is also shown in 
Fig. O Normalizing this to the data, the total cross sec- 
tion was determined to be 3.99 ± 0.61(stat) ± 0.58(sys) 
± 0.40(abs) fib. The quoted systematic error of 14% was 
estimated by setting the measured cross sections all to 
their upper systematic error limits or all to their lower 
systematic error limits and noting how the cross section 
changed. The variation in the total cross section ex- 
tracted if we use the same procedure with different PDF 
choices and models was estimated to be small ( 3%). 

A comparison is made in Fig. 0Ji of the present (pt) 
value with values from previous experiments 0. There 
are no theoretical predictions that we can compare with 
the (pt) measurements. The total J/4> cross section 
determined in this analysis is shown in Fig. 0Jd, along 
with cross sections determined by lower energy experi- 
ments and predictions from the COM |^ using two 
different PDFs. The y/s dependence of the cross section 
is sensitive to the factorization scale Q, since the shape 
of the PDFs depend on Q. The values of Q (3.1 GeV 
for GRV98NLO and 2.3 GeV for MRST2001NLO) were 
chosen to give good agreement with the data. The to- 
tal cross section normalization was obtained using color 
octet matrix elements from but has large theo- 

retical uncertainties associated with the charm quark 
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FIG. 3: The central rapidity point is from J/ip — > e'^e~ , the 
others from J/i/i /i^ jj.^ . The brackets represent systematic 
uncertainties. All curves have their overall normalization fit- 
ted to the data. The Pythia shape was used to determine the 
cross section. There is an overall 10% absolute normalization 
error not shown. 
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FIG. 4: (a) The present J/tp mean pr value compared with 
previous measurements at lower energy. The linear fit param- 
eters are p = 0.53, q = 0.19. (b) The present J/ip total cross 
section compared with previous measurements at other values 
of ^/s. The curves are discussed in the text. 

mass and the renormalization scale. The renormaliza- 
tion scale was taken to be equal to the quark mass Mc, 
and their values (1.48 GeV for GRV98NLO and 1.55 GeV 
for MRST2001NLO) were chosen to give good agreement 
with the data. The CEM is also able to describe the to- 
tal cross section data 6J. All measurements and models 
include feed-down from the Xc E^nd the V'' to the J/ijj. 
We estimate (l^i] that B decay feed-down contributes less 
than 4% to the J/ tp total cross section at y/s = 200 GeV. 

In summary, we have presented the first pp J/^ + X 
measurements from RHIC, obtained at ^/s = 200 GeV. 
The rapidity distributions, px distributions, (pt) and 
total cross sections have been presented and compared 
with available model calculations. The transverse mo- 
mentum distributions above 2 GeV/c are reasonably well 
described by the COM. With the present statistical pre- 
cision, our rapidity distribution shape is consistent with 
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most of the available PDFs. COM calculations are able 
to reproduce the ^/s dependence of the cross section us- 
ing color octet matrix elements found in the literature, 
with a reasonable choice of QCD parameters. 

RHIC is expected to have proton-proton runs with en- 
hanced luminosity at y/s — 200 and 500 GeV in the near 
future. The increased luminosity will improve the statis- 
tical precision and pr reach of the PHENIX data, and will 
ultimately make it possible to measure the J/ ip polariza- 
tion, which has been an important test for models 
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